Converse Magnetoelectric Effects in FesOzt/BaTiOs Multiferroic Hybrids 
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The quantitative understanding of converse magnetoelectric effects, i.e., the variation of the mag- 
netization as a function of an applied electric field, in extrinsic multiferroic hybrids is a key pre- 
requisite for the development of future spintronic devices. We present a detailed study of the 
strain-mediated converse magnetoelectric effect in ferrimagnetic FesCU thin films on ferroelectric 
BaTiOs substrates at room temperature. The experimental results are in excellent agreement with 
numerical simulation based on a two-region model. This demonstrates that the electric field induced 
changes of the magnetic state in the FesCU thin film can be well described by the presence of two 
different ferroelastic domains in the BaTiOs substrate, resulting in two differently strained regions 
in the Fes04 film with different magnetic properties. The two-region model allows to predict the 
converse magnetoelectric effects in multiferroic hybrid structures consisting of ferromagnetic thin 
films on ferroelastic substrates. 

PACS numbers: 75.70.Cn 75.80. +q 75.85. +t 81.15.Fg, 85.75.-d 85.80.Jm 
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I. INTRODUCTION 

The manipulation of the magnetization by electric 
fields, which is referred to as converse magnetoelectric 
effect, is one of the most challenging tasks in today's spin- 
tronic devices. 1 3 The presence of such a cross-coupling 
between ordered magnetic and dielectric states in so- 
called multiferroic systems allows for the development of 
a new family of devices, where the spin degree of freedom 
is controlled by electric instead of magnetic fields. 4 10 It 
has been shown that a robust electric field control of mag- 
netism can be realized by using extrinsic multiferroic hy- 
brid structures consisting of ferroelectric and ferromag- 
netic materials. 11 15 In these systems, the extrinsic mag- 
netoelectric effects rely on either electric field effects 16 
using carrier mediated ferromagnets, 17 or employing ex- 
change coupling effects between antiferromagnetic, fer- 
roelectric and ferromagnetic compounds. 18 ' 19 Moreover, 
from the beginning the elastic coupling at the interface 
between ferromagnetic and ferroelectric materials has 
been exploited to realize a robust electric field control 
of the magnetization. 20 22 In this case, the electric and 
magnetic degrees of freedom are elastically coupled at the 
interface between the two ferroic constituents. More pre- 
cisely, strain mediated converse magnetoelectric effects 
rely on an electric field induced strain in the ferroelectric 
via the converse piezoelectric effect, or via ferroelastic 
reorientations. This strain is then transferred into the 
ferromagnetic thin film clamped onto the ferroelectric, 
where converse magnetoelastic effects modify the mag- 
netization. 

Indirect strain-mediated magnetoelectric effects 
can be realized using particulate composites, 23 lam- 
inated composites, 24 nanostructured composites, 25 
and horizontal multiferroic hybrid structures. In 
the latter case, the discovery of giant sharp and 
persistent converse magnetoelectric effects in ferromag- 



netic thin films epitaxially grown on BaTiOs (BTO) 
crystals 22 triggered an increasing research effort in this 
field, using different ferromagnetic material systems, 
such as magnetic ferrites (Fe 3 4 , 26 31 CoFe 2 4 , 32 ' 33 
and NiFe204 34 ), magnetic perovskites and double- 
perovskites (Lai-^Sr^MnOs, 35 Lai-^Ca^MnOs, 36,37 
Pri^Ca^MnOa, 38 and Sr 2 CrRe0 6 39 41 ), and 3d- 
magnets (Fe, 42 45 Ni, 5 ' 46 49 Co, 50 and CoFe 51 ). In 
particular, strong electric field-induced changes of the 
magnetic state were reported in Fe304-based extrinsic 
multiferroic hybrids. 26,52 54 Since Fe304 has a ferrimag- 
netic ground state with a high Curie temperature of 
about 850 K, 55 a pronounced magnetostrictive effect, 56 
and a high spin polarization at the Fermi level, 57 
Fe304 is a promising candidate for possible future 
magnetoelectric devices. Moreover, a sizable sponta- 
neous ferroelectric polarization was reported in Fe304 
films 58,59 which is caused by a non-centrosymmetric 
charge order. 60,61 Therefore, Fe304 seems to be not 
only the oldest magnetic material but also the first 
magnetoelectric multiferroic known to mankind. 62 

To quantitatively understand the strain-mediated 
magnetoelectric effects in BTO-based multiferroic hybrid 
structures, as a first step, we recently investigated the 
magnetization changes as a function of temperature in 
FesoCoso/BTO and Ni/BTO hybrids caused by magne- 
toelastic effects. 63 We demonstrated that the manipula- 
tion of the magnetization in BTO-based multiferroic hy- 
brid structures can be theoretically well described by con- 
trolling the ferromagnetic and ferroelastic domain config- 
uration in the multiferroic hybrid. Here, as a second step, 
we show that by knowing the ferroelastic domain config- 
uration in the BTO substrate as a function of the ap- 
plied electric field, the strain-mediated converse magne- 
toelectric effect in epitaxially grown Fes04/BTO hybrid 
structures can be theoretically simulated using a simple 
model, consisting of two magnetic regions in the Fes04 




FIG. 1. (color online) Intensity evolution of the RHEED (00) reflection monitored during the deposition of a 64 nm thick FesCU 
film on a 0.5 mm thick BTO substrate. The intensity was integrated within the red rectangle marked in the RHEED patterns, 
which are shown in the insets. The time stamps at which the RHEED patterns were recorded are marked by (I)-(IV). The 
growth process can be subdivided into three main steps: highly disordered growth (A) (0 s < t < 25 s) — three dimensional 
growth mode (B) (300 s < t < 3000 s) — layer-by layer growth mode (C) (3000 s < t < 3950 s). 



thin film. This demonstrates that by controlling the fer- 
roelastic domain state in the BTO substrate large and 
robust manipulations of the magnetization as a function 
of the applied electric field are feasible at room temper- 
ature in BTO-based multiferroic hybrid structures. 

The article is organized as follows: In Sec. II, we 
discuss the growth process and the structural as well 
as magnetic properties of FesCU thin films epitaxially 
grown on BTO substrates. Since X-ray diffraction mea- 
surements using synchrotron radiation carried out on 
these Fes04/BTO hybrid structures reveal two differ- 
ently strained regions in the Fe304 thin film at low elec- 
tric fields, a two-region model is introduced in Sec. III. In 
this model, the room-temperature magnetization behav- 
ior of Fes04/BTO hybrids as a function of the applied 
electric field can be explained by the magnetic behav- 
ior of parts of the Fe304 thin film elastically clamped 
onto two different ferroelastic domains of the BTO sub- 
strate. Thus, to simulate the strain-mediated converse 
magnetoelectric effect, at first the volume fraction of the 
ferroelastic domains in the BTO substrate is determined 
(cf. Sec. IV), and then the magnetic state of the Fes04 
thin film on top of these domains is calculated using mag- 
netoelastic theory (cf. Sec. V). In Sec. VI, we show that 
there is an excellent agreement between the experimen- 
tal results and the theoretical simulations based on the 
two-region model. This demonstrates that the experi- 
mentally obtained room-temperature converse magneto- 
electric effects in BTO-based multiferroic hybrids are well 
described within this approach for magnetic fields larger 
than the magnetic saturation field. 



II. THIN FILM GROWTH 

The Fes04 thin films with thicknesses between 40 nm 
and 70 nm were epitaxially grown on 0.5 mm thick (001)- 



oriented BTO substrates by laser molecular-beam epi- 
taxy ( laser- MB E) 64 monitored by an in-situ high pres- 
sure RHEED system. 65 ' 66 The energy density of the KrF 
excimer laser (A = 248 nm) at the target was set to 
3.1J/cm 2 and the laser repetition rate was 2 Hz. The 
deposition was carried out in an atmosphere of pure ar- 
gon with a pressure of 1.6 x 10 -3 mbar at a temperature 
of 593 K, i.e., in the cubic phase of the BTO substrate. 

The growth process of the Fe304 thin films is illus- 
trated in Fig. 1 on the basis of the intensity evolution 
of the RHEED (00) spot, which was recorded during the 
growth of a 64 nm thick Fe304 film. Before starting the 
growth process, the RHEED pattern of the BTO sub- 
strate reveals three spots, which belong to the first Laue 
circle and can be indexed with (11), (00), and (11) [cf. in- 
set (I) of Fig. 1]. By starting the deposition at t = s, 
the high intensity of the (00) spot drastically decreases 
[cf. Fig. 1 (A)]. After 50 pulses (t « 25s), the deposi- 
tion was stopped [cf. Fig. 1 (II)]. The RHEED pattern 
recorded at this moment does not show any reflections 
[cf. inset (II) of Fig. 1]. This indicates a highly dis- 
ordered and nearly amorphous growth within the first 
monolayers of Fes04. Further deposition leads to a three 
dimensional crystalline growth as indicated by a checker- 
board like RHEED pattern [cf. inset (III) of Fig. 1] . Upon 
further continuing the deposition, the three-dimensional 
growth mode turns over into a two-dimensional one with 
increasing RHEED intensity and emerging intensity os- 
cillations [cf. Fig. 1 (B)]. In the last stage of the growth 
process (t > 3000 s), a real two-dimensional layer-by- 
layer growth mode can be obtained [cf. Fig. 1 (C)]. At 
this stage, a single unit cell of Fe304, which is manifested 
by four RHEED oscillations, 67 was deposited followed by 
a growth interruption of 30 s, allowing for the relaxation 
of the film surface. The RHEED pattern recorded at 
the end of the growth process [cf. inset (IV) of Fig. 1] 
does not show any indication of transmission like spots, 
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FIG. 2. (color online) Magnetic hysteresis of a 45 nm thick 
Fe3 04 film grown on a (OOl)-oriented BTO substrate mea- 
sured at 300 K with the magnetic field applied in the film 
plane. The magnetic field dependence of the magnetiza- 
tion can be fitted by a (1 — b/y/H) behavior using b = 
(0.045±0.005) VT (red solid line), (a) The M(H) loop reveals 
a coercive field of 35 mT and a saturation magnetization of 
M s = 438kA/m at [iR = 7T. (b) The normalized remanent 
magnetization as a function of temperature measured after 
field cooling the sample with iiqH — 7T discloses a Verwey 
transition at 121 K. 



suggesting a smooth surface of the Fe3C>4 thin film. Us- 
ing atomic force microscopy as well as X-ray reflectom- 
etry, a surface roughness (rms value) ranging between 
0.8 nm and 2.5 nm was found. Therefore, as obvious from 
Fig. 1, Fes04/BTO hybrid structures can be fabricated in 
a true layer- by- layer growth mode. However, the highly 
disordered growth at the beginning of the deposition of 
Fes04, which might be explained by the lattice mismatch 
of Fe304 and BTO of 4.7% at the growth temperature, 
results in a large mosaic spread. X-ray diffraction mea- 
surements (not shown here) reveal no secondary phases 
as well as a full width at half maximum of the rocking 
curves around the FesC^ (004) reflection of about 0.6°. 
This is more than one order of magnitude larger than 
observed for Fe304 thin films grown on lattice-matched 
MgO substrates. 68 

The magnetic properties of the hybrid structures were 
studied via superconducting quantum interference device 
(SQUID) magnetometry in the temperature range be- 
tween 2 K and 300 K at magnetic fields of up to 7 T ap- 
plied in the film plane. Figure 2 shows the magnetic 
behavior of a 45 nm thick Fe304 film grown on a BTO 
substrate. The magnetic hysteresis loop M(H) measured 
at 300 K reveals a coercive field of 35 mT [cf. inset (a) of 
Fig. 2], which is close to the value reported for Fes04 thin 
films on lattice-matched MgO substrates. 69 Furthermore, 
a high saturation magnetization of M s = 438kA/m at 
fioH = 7T is obtained, which corresponds to 3.1 /ie/f-u.. 
The difference AM S between M s and the bulk value 
M s bulk = 4.1/iB/f-u. is generally explained by the pres- 



ence of anti-phase boundaries (APB). APBs exhibit a 
strong antiferromagnetic coupling, which results in a 
(1 — b/y/H) behavior of the magnetization, where the 
parameter b depends on the density of APBs. 70 Using 
b = (0.045 ± 0.005) a/T, the magnetic hysteresis can be 
nicely fitted [cf. red line in Fig. 2]. This value is much 
lower than values reported for Fe304 thin films grown on 
(OOl)-oriented MgO or A1 2 3 substrates. 70 ' 71 This might 
be explained by the disordered growth within the initial 
phase of the deposition process, leading to a disconnec- 
tion of the Fes04 thin film with the surface structure of 
the BTO substrate. Thus, the deviation in saturation 
magnetization AM S might not be explained by APBs 
alone. Oxygen non-stoichiometry Fe 3 ( 1 _5)04 with S ^ 
is another possibility, which affects M s . 72 The degree of 
oxygen non-stoichiometry S can be investigated by mea- 
suring the Verwey transition Ty, which is a hallmark of 
the quality of Fe304 thin films. The inset (b) of Fig. 2 
shows the remanent magnetization as a function of tem- 
perature measured after field cooling the sample in an 
external magnetic field of fi^H — 7T. The Verwey tran- 
sition is clearly observable at 121 K, which is identical 
to the bulk value. Thus, our Fes04 thin films grown 
on BTO substrates have a low amount of oxygen vacan- 
cies and the correct iron oxide phase. Other iron oxides 
such as maghemite can be safely excluded. Note that 
a small reduction of the phase transition temperature 
was observed in other samples, which can be translated 
into an oxygen non-stoichiometry of 5 = 0.0021 assum- 
ing a linear dependence of Ty on S. 73 Altogether, the 
measured reduction of the saturation magnetization with 
respect to the bulk value might be explained by a com- 
bination of effects caused by APBs and a slight oxygen 
non-stoichiometry in the Fe304 thin films. 



III. TWO REGION MODEL 

One of the intriguing properties of multiferroic hybrid 
structures consisting of ferromagnetic and ferroelectric 
materials is the possibility to manipulate the magnetiza- 
tion M by an electric field E. To investigate this converse 
magnetoelectric effect in the Fes04/BTO hybrid struc- 
tures, an Au bottom electrode was sputtered on the back- 
side of each BTO substrate, which enables us to apply 
an electric field E z across the BTO substrate along the z 
direction using Fe304 as top electrode. The Fes04/BTO 
hybrids are further heated to 450 K well above the Curie 
temperature of BTO and slowly cooled down to room 
temperature while applying an electric field of 400 kV/m. 

As an example, Fig. 3 shows the modification of 
the magnetization as a function of the applied electric 
field, M{E), for different values of the applied mag- 
netic field. The curves were measured at room tem- 
perature, i.e., in the tetragonal phase of BTO, using a 
Fe304(38nm)/BTO multiferroic hybrid. Here, the pro- 
jection of the magnetization M on the magnetic field 
direction, which was oriented along the y direction, is 
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FIG. 3. (color online) Converse magnetoelectric effects mea- 
sured at 300 K for a Fe 3 4 (38nm)/BTO multiferroic hybrid. 
The magnetic field was applied in the film plane along the 
y direction. A miscut BTO substrate was used for the ex- 
periment with a misalignment of 1° (7°) of the [001] ([100]) 
direction of the ferroelastic unit cell with respect to the z (x) 
direction describing the surface of the BTO substrate. The 
variation of the ferroelastic domain configuration as a func- 
tion of the electric field is illustrated in the bottom part of 
the figure. At high positive and negative fields, a ferroelastic 
single c domain state is expected, while at low electric fields 
a c/d2 multi-domain state is present in the BTO crystal. The 
blue arrows indicate the orientation of the ferroelectric polar- 
ization and the blue regions denote ferroelastic domain walls. 
Polarization rotation induced by the electric field is neglected 
for simplicity. 



plotted as a function of the electric field strength E z ap- 
plied across the BTO substrate. In analogy to Ni/BTO 
hybrid structures discussed in Ref. 46, a butterfly shape 
of the M(i£)-loops is observable. On the basis of the 
measured magnetization changes AM, a pseudo-linear 
magnetoelectric coupling coefficient can be estimated us- 
ing a 2 3 = <^fm/(^fm + cIfe)^oAM/AE, which results 
in Qf23 = 3 x 10 _15 s/m. This value is lower than the 
magnetoelectric constant of the prototype intrinsic mag- 
netoelectric material C^Os (a = 4.13 x 10 _12 s/m), 74 
demonstrating that the magnitude of the converse linear 
magnetoelectric effect in BTO-based hybrid structures is 
comparable or even less than in intrinsic magnetoelec- 
tric materials due to the unfavorable ratio between the 
thickness of the ferromagnetic thin film (g^fm) and the 
thickness of the BTO substrate (c?fe)- However, the 
observed magnetoelectric effect is robust at room tem- 
perature. Only very recently, large room-temperature 
magnetoelectric effects were reported in magnetoelectric 
hexaferrites. 75,76 This was the first observation of signif- 
icant magnetoelectric effects in intrinsic magnetoelectric 
materials at room temperature. 
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FIG. 4. (color online) X-ray diffraction around the asym- 
metric BTO (204) reflection of a Fe 3 4 (38nm)/BTO hybrid 
recorded in the (a) qHoo— gooz,-plane and (b) qoKo— QooL-plane 
at OkV/m using a miscut BTO crystal (yellow/red: high in- 
tensity, light blue/dark blue: low intensity). Only reflections 
caused by ferroelastic c domains and tilted a^ domains are 
visible. The expected positions of the reflections caused by 
diffraction from tilted ferroelastic a\ domains are marked by 
dashed circles. 



The variation of the magnetization as a function of an 
applied electric field in BTO based multiferroic hybrids 
at 300 K can be explained by taking into account the 
ferroelastic domains of BTO. 46 In the tetragonal phase, 
BTO exhibits three types of ferroelastic domains with 
the polarization vector pointing either along the out-of- 
plane direction (c domains) or one of the two in-plane 
directions [a\ domains and a 2 domains). The tetrago- 
nal unit cell of the a domains is oriented with its longer 
axis parallel to the BTO surface. It induces tensile strain 
in the overlying ferromagnetic thin film, which modifies 
the magnetic anisotropy due to converse magnetoelastic 
effects. Thus, to quantitatively describe converse mag- 
netoelectric effects in BTO-based multiferroic hybrids, 
the volume fraction of the different ferroelastic domains 
has to be known. The requirements are relaxed by using 
miscut BTO crystals. In this case only a single type of a 
domains is formed. 

Figure 4 shows x-ray diffraction measurements around 
the asymmetric BTO (204) reflection. In the experiment 
a Fe 3 4 (38nm)/BTO hybrid with a BTO substrate is 
used, which has a misalignment of about 1° (7°) of the 
[001] ([100]) direction of the ferroelastic unit cells with 
respect to the z (x) direction describing the surface of 
the BTO crystal (cf. inset of Fig. 3). The mesh scan 
around the BTO (204) reflection reveals only two reflec- 
tions. One is caused by scattering from ferroelastic c do- 
mains at qooL = 4.0 rlu and the other by tilted ferroelastic 
d2 domains. The second type of a domains (a\ domains) 
could not be detected, neither in the qnoo — gooL-plane 
[cf. Fig. 4(a)] nor in the q xo ~ gooL-plane [cf. Fig. 4(b)]. 
Therefore, only a single type of a domains is present in 
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FIG. 5. (color online) (a) X-ray diffraction around the BTO 
(002) reflection recorded at an electric field of 400kV/m 
(black full symbols) and subsequently of -28 kV/m (red full 
symbols) applied across a FesC^/BTO hybrid structure. The 
measurements were performed at the BM28 beamline at the 
European Synchrotron Facility (ESRF). Upon increasing the 
electric field from -28kV/m to 400kV/m again (black open 
symbols), the initial measurement at 400kV/m is reproduced. 
The effect of the strain imposed on the FesC^ film is shown 
in more detail in (b) and (c). In the case of 400kV/m, the 
(004) reflection of the FesCU thin film (black symbols) can be 
well fitted using a single peak profile (black solid line). At 
kV/m, two peak profiles (dashed lines) are needed to fit the 
measured reflection (red solid line) . The gray lines display the 
difference between the experimental data (symbols) and the 
fits (solid lines), and the blue arrows mark the gooz, positions 
of the reflections. 



the used miscut BTO crystals at zero electric field. In 
this case the evolution of the ferroelastic domains as a 
function of the applied electric field E z can be described 
as follows: At high positive and negative electric fields, 
a ferroelastic single c domain state is present, which 
changes into a ferroelastic multi-domain state around 
E z — OkV/m. This ferroelastic multi-domain state con- 
sists of ferroelastic c and tilted <22 domains as well as cc 
and ac domain walls (cf. bottom part of Fig. 3). 

The impact of the ferroelastic domain state on the 



overlying ferromagnetic thin film due to the imposed 
strain can be investigated by performing x-ray diffraction 
measurements on crystalline FeaO^BTO hybrid struc- 
tures. As an example, Fig. 5 shows L-scans around 
the BTO (002) reflection recorded at different electric 
field strengths E z applied across the FesO^BTO hy- 
brid structure. The measurements shown in Fig. 5 
were carried out at the BM28 beamline at the Euro- 
pean Synchrotron Radiation Facility (ESRF) using syn- 
chrotron radiation with an energy of 7.1255 keV. At 
E z — 400kV/m, the BTO substrate yields only one 
reflection at ^ool = 2.0 rlu [cf. full black symbols in 
Fig. 5(a)]. This confirms that a single ferroelastic c do- 
main state exists at this electric field strength. By reduc- 
ing the electric field to values close to the coercive field 
of BTO (E z = — 28kV/m), a second reflection around 
Qool = 2.03 rlu emerges, which is caused by ferroelastic 
d2 domains formed in the BTO substrate [cf. full red sym- 
bols in Fig. 5(a)]. By increasing the electric field to the 
starting value of E z = 400 kV/m, the first measurement 
is reproduced, i.e., a single ferroelastic c domain state is 
restored [cf. open black symbols in Fig. 5(a)]. The im- 
pact of the strain imposed on the Fes04 thin film was 
investigated in detail around the Fes04 (004) reflection 
[Figs. 5(b) and (c)]. The peak profile of the Fe 3 4 (004) 
reflection measured while applying an electric field of 
E z = 400kV/m can be well fitted using a single Pear- 
son VII peak function 77 [cf. solid line in Fig. 5(b)]. This 
demonstrates that the Fes04 thin film is homogeneously 
strained at this electric field strength. An out-of-plane 
lattice constant of c = (0.8492 ± 0.0005) nm can be de- 
rived from the peak position of the Fes04 (004) reflection. 
Furthermore, x-ray diffraction measurements carried out 
around the Fe304 (404) reflection (not shown here) reveal 
an in-plane lattice constant of a = (0.8296 ± 0.0008) nm. 
Comparing these values to the bulk lattice parameters in- 
dicates that our Fes04 thin films exhibit a tensile out-of- 
plane strain of +1.2% and a compressive in-plane strain 
of —1.2%. This is in contrast to x-ray diffraction mea- 
surements published recently. 26 ' 29 However, those mea- 
surements were performed without any electric field ap- 
plied to the samples, i.e., in a ferroelastic multi-domain 
state of the BTO substrate. The derived values of the in- 
plane and out-of-plane strain yield a Poisson ratio of 0.33, 
which is in good agreement with Ref. 78. In contrast to 
the x-ray measurements at E z = 400 kV/m, the intensity 
of the Fe 3 4 (004) reflection at E z = OkV/m is strongly 
reduced and a second peak occurs close to the first one 
[Fig. 5(c)]. Two different peak functions are needed to 
reproduce this situation [cf. dashed lines in Fig. 5(c)]. 
While the position of the first one is almost identical to 
the peak position measured at E z = 400 kV/m, the sec- 
ond peak is located at a higher gooL-value [cf. blue arrows 
in Figs. 5(b) and (c)]. This clearly demonstrates that two 
differently strained regions are present in the Fes04 thin 
film at E z = OkV/m. The one region can be attributed 
to those parts of the Fe304 thin film clamped to ferroe- 
lastic c domains of the BTO crystal, while the other to 



those regions of the FesCU thin film positioned on top of 
ferroelastic ct2 domains. The <22 domains induce a ten- 
sile strain along the y direction in the overlaying Fe304 
thin film, thereby causing a compressive strain along the 
out-of-plane direction due to elasticity. This shifts the 
Fe304 (004) reflection to higher #ool values. From the 
difference in position of the two peaks visible in Fig. 5(c), 
the induced compressive out-of-plane strain of the FeaC^ 
thin film due to the transformation of ferroelastic c do- 
mains into <22 domains in the BTO crystal can be cal- 
culated to —(0.53 ± 0.03)%. By employing the Poisson 
ratio of 0.33, the induced in-plane strains e^? = 0% and 
e l yy = (1.07 ±0.06)% are derived. These values are in ex- 
cellent agreement with the theoretical values e^J° = 0% 
and Cyy° = 1.05% caused by domain reorientations from 
ferroelastic c to a^ domains in the BTO crystal. This 
provides clear evidence for a perfect strain transmission 
from the BTO substrate into the Fes04 thin film. 

Figure 5 indicates that the large tensile strain is im- 
posed on the Fe304 thin film by the formation of ferroe- 
lastic <22 domains in the BTO substrate. Due to mag- 
netoelastic effects this strain should strongly affect the 
magnetization M a of those parts of the Fes04 film elas- 
tically clamped onto these domains. In contrast, since 
the (004) reflection of the Fe304 thin film positioned on 
top of ferroelastic c domains stays nearly unaffected as 
a function of the applied electric field, the magnetiza- 
tion M c of these regions is expected to depend only very 
weakly on the applied electric field. Thus, the converse 
magnetoelectric effects in FeaO^BTO hybrid structures 
at room temperature can be described in first order by a 
two-region model 
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FIG. 6. (color online) Determination of the volume fraction of 
ferroelastic c domains x c (E) using x-ray diffraction, (a), (b) 
Mesh scans around the BTO (002) reflection recorded while 
applying an electric field of 400 kV/m and kV/m across the 
Fe3 04(38nm)/BTO hybrid. To obtain x c , the ratio of the 
integrated intensities for gooL < 2.01 rlu and gooz, > 2.01 rlu 
is calculated, (c) x c as a function of the applied electric field. 



Here, x c (E) denotes the volume fraction of the ferroelas- 
tic c domains in the BTO substrate as a function of the 
applied electric field, and (1 — x c (E)) the remaining vol- 
ume fraction of the emerging ferroelastic a^ domains. In 
this model, the magnetizations M a and M c are consid- 
ered independent of the electric field. Thus, Eq. (1) only 
considers ferroelastic domain reorientations and does not 
take into account linear converse piezoelectric effects as 
well as effects caused by polarization rotations. 79 This is 
reasonable, since the inverse linear piezoelectric effect of 
BTO leads to a maximum in-plane strain, which is three 
orders of magnitude smaller than the strain caused by 
ferroelastic domain reorientation. 46 Furthermore, effects 
due to ferroelastic domain walls are also not included in 
the model. To calculate the converse magnetoelectric ef- 
fects using Eq. (1), the volume fraction x c (E) as well as 
the magnetization values M c and M a have to be deter- 
mined. 



IV. VOLUME FRACTION OF FERROELASTIC 
DOMAINS 

The volume fraction of the ferroelastic c domains 
x c (E) is determined by means of x-ray diffraction 
measurements. 80 Since the ferroelastic c and a domains 
exhibit different unit cell dimensions as well as a change 
of the polarization direction, both the Bragg angles and 
the intensities of the respective reflections are affected. 81 
For the experimental determination of x ci x-ray diffrac- 
tion measurements around the BTO (002) reflection were 
carried out for different values of the electric fields ap- 
plied across the Fe 3 04/BTO hybrid structure. 82 Fig- 
ures 6(a) and (b) show two reciprocal space maps 
recorded with E z = 400kV/m and E z = OkV/m. 
For E z = 400kV/m, only a single reflection around 
Qool — 2.0 rlu is visible, which is a marker for ferroe- 
lastic c domains. This demonstrates that at this elec- 
tric field strength a purely ferroelastic c domain state is 



present in the BTO substrate. By reducing the electric 
field strength, tilted a 2 domains are formed continuously, 
until the electric coercive field is reached. Scattering 
from these <22 domains causes a finite intensity around 
q 00L = 2.02 rlu [cf. Fig. 5(b)]. To determine x c (E), the 
detected intensities were integrated for #ool < 2.01 rlu 
(J c ) and <2ool > 2.01 rlu (J a ). Using these intensities, 
x c (E) can be estimated as x c (E) = I c /(I a + I c ). The re- 
sulting volume fraction as a function of the applied elec- 
tric field is shown in Fig. 6(c). By comparing x c (E) with 
the M(^)-loops depicted in Fig. 3, the pronounced in- 
terdependence of the magnetization in the ferromagnetic 
thin film and the volume fraction of the ferroelastic c do- 
mains in the BTO substrate becomes obvious. This pro- 
vides further evidence for the validity of the two-region 
model described by Eq. (1). We note, however, that for 
the calculation of M(E), the ferroelastic domain configu- 
ration at the surface of the BTO is essential, since these 
domains are elastically coupled to the overlaying ferro- 
magnetic thin film. As X-rays penetrate deep into the 
BTO substrate, the obtained volume fraction x c (E) is an 
averaged value. By using different angles of the incom- 
ing x-ray beam, i.e., varying the penetration depth of the 
X-rays, an uncertainty of less than 3% for x c (E) can be 
estimated. As obvious from Fig. 6(c), this uncertainty is 
not larger than the already displayed error bars, describ- 
ing the experimental uncertainty at a constant angle of 
the incoming X-rays. 



magnetization M s and assume M s to be homogeneous 
throughout the Fe304 thin film. In this approximation, 
the magnetization values M c and M a are described by 
M c = M s m c and M a = M s m a , respectively. Using mag- 
netoelastic theory, the unit vectors m c and m a can be de- 
rived by minimizing the free enthalpy density g FM of the 
Fe304 thin film, which is based on the magnetic energy 
density u FM given by (see Ref. 63 for more details) 
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Here, ^™( m * s the magnetic anisotropy contribution, 
^™(^/e) the purely elastic energy density, ^™ gel the 
magnetoelastic energy density, and rrii with i = 1, 2, 3 are 
the components of the unit vectors m a>c , i.e., the direc- 
tional cosines. The strain components rjk (k — 1, . . . , 6) 
are given in matrix notation with 771 = en, 772 — ^22, 
m = ^33, m = 2e 23 , rj 5 = 2e 3i, and rj 6 = 2ei 2 . 84 To 
account for the shape anisotropy in ferromagnetic thin 
films with finite dimensions, an additional contribution 
«demag(^3) = (Mo/2)M s 2 m§ is added to Eq. (2). 

In the framework of a magnetic single domain model, 
u FM is only a function of the magnetization direction rrii 
and the strain components rjk. The components rjk with 
k > 4 are zero, since no shear strains are expected in 
the tetragonal phase of BTO at room temperature. As 
discussed in Section III, the strain state of those parts of 
the Fes04 thin film elastically coupled to ferroelastic c 
domains (r] c ) and a domains (rj a ) of the BTO substrate 



V. MAGNETOELASTIC EFFECTS 

Having determined the ferroelastic domain state in 
the BTO substrate as a function of the applied elec- 
tric field strength, in the two-region model approach 
[cf. Eq. (1)] the strain- mediated converse magnetoelec- 
tric effects M(E) in the Fe 3 4 /BTO hybrid depend only 
on the magnetization values M c and M a of the regions 
elastically clamped to ferroelastic c and <22 domains, re- 
spectively. Therefore, the remaining step is to determine 
the magnetization values M c and M a . 

In general, elastic strain can modify the magnetization 
by two main effects. First, the direction of the magnetiza- 
tion can be altered by a strain-induced change of the mag- 
netic anisotropy, since the magnetoelastic energy can be 
regarded as an additional uniaxial magnetic anisotropy 
in a phenomenological free energy approach. 5 Second, 
due to strain-induced changes of the bond length and 
bond angle, the exchange interaction can be modified. 
This may cause variations of the magnetic ordering tem- 
perature and the magnetization. Even a strain-induced 
transition from a magnetically ordered/electrically con- 
ducting to a nonmagnetic/insulating state is possible by 
changing the spin state of the magnetic ions. 83 Since 
previous experiments showed that changes of the mag- 
netic anisotropy are the driving force behind the converse 
magnetoelectric effects in FeaO^BTO hybrids, 63 we ne- 
glect any strain-induced modification of the saturation 
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Assuming a cubic symmetry of the Fes04 thin 
film, the magnetic anisotropy ^™( m can be ex ~ 
pressed as u™(rrii) = K c (m\m^ + m^m\ + m^rri^), 
with the first order cubic anisotropy constant K c . 
The second term in Eq. (2) represents the purely 
elastic energy ^™fe) = |cn (rjf + r\\ + rjf) + 

C12 (^71772 + ^2^3 + ^1^3) + 2^44 (^ + ^5 + Ve)' Here > C ij 

are the elastic stiffness constants of the cubic Fe304 thin 
film, which are regarded as material constants (en = 
27.2 x 10 10 N/m 2 , c 12 = 17.8 x 10 10 N/m 2 , and c 44 = 
6.1 x 10 10 N/cm 2 ). 78 The interaction between the elastic 
and the magnetic anisotropy energies is described by the 
first order magnetoelastic energy density ^™ g ei : 
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FIG. 7. (color online) Magnetization M plotted versus the an- 
gle <j> denning the in-plane orientation of the applied magnetic 
field. = 0° corresponds to the x direction. The data are 
taken at an electric field of E z — 400 kV/m and an external 
magnetic field of 25 mT (black symbols) and 20 mT (red sym- 
bols). The M(0)-curves can be fitted using a cubic anisotropy 
field of K c /M s = (— 19 =b 1) mT (solid lines), assuming a mis- 
alignment of 4° between the sample normal and the rotation 



with the magnetoelastic coupling coefficients B\ and B 2 . 
For bulk magnetostrictive samples, B\ and B2 can be ex- 
pressed as a function of the magnetostrictive strains A100 
and Am: B 1 = -§Aioo (en - c 12 ) and B 2 = -3A m c 4 4. 
We note that the magnetoelastic coupling coefficients 
Bi in ferromagnetic thin films may deviate from the 
bulk values due to surface effects and/or the influence 
of strain. 85 ' 86 In our analysis, we neglect these effects 
and use the bulk values A100 = —19.5 x 10 -6 and Am = 
+77.6 x 10 -6 for the calculations. 87 A proportionality 
factor x is introduced to account for any deviation in the 
magnetoelastic coupling from bulk-like behavior. 

To determine the cubic anisotropy field K C /M S1 
SQUID magnetometry measurements were carried out 
for different in-plane orientations of the applied mag- 
netic field. To ensure a well defined ferroelastic single 
c domain state in the BTO substrate, an electric field of 
E z = 400 kV/m was applied across the hybrid structure. 
The angular dependence of the magnetization projection 
M((j))/M((j) = 0°) along the external magnetic field H is 
shown in Fig. 7 for fi H = 20 mT and ji Q H = 25 mT. Ob- 
viously, a four- fold symmetry is visible, which is a clear 
sign of a cubic magnetic anisotropy in the Fe304 thin 
film. The angular dependence M{(j))/M{(j) = 0°) mea- 
sured at constant strain state r] c can be simulated by 
means of Eq. (2). The best fit between experiment and 
simulation was obtained using a cubic anisotropy field 
of K c /M s = (-19 ± l)mT (cf. solid lines in Fig. 7). A 
misalignment of 4° between the sample normal and the 
rotation axis was assumed to account for the different 
magnetization values measured at <\> ~ —80°, (j) ~ 10°, 
and <j>^ 100° (cf. Fig. 7). 
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FIG. 8. (color online) Measurement (symbols) and sim- 
ulation (lines) of the converse magnetoelectric effect in a 
Fe3 04(38 nm)/BTO multiferroic hybrid for different magnetic 
field strengths. The magnetic field was aligned in the film 
plane along the y direction of the BTO crystal and the tem- 
perature was set to 300 K. For better visibility, the curves are 
vertically shifted with respect to each other. The best fit be- 
tween experiment and simulation was obtained by assuming a 
magnetoelastic coupling efficiency of 55% (x = 0.55). A mis- 
alignment of the sample surface with respect to the magnetic 
field of 5° was assumed for the calculation. 



Using the cubic anisotropy field and the elastic param- 
eters determined above, the FesO^BTO hybrid struc- 
ture is well described except for the remaining unknown 
parameter %, which characterizes the magnetoelastic cou- 
pling strength. This parameter is the only fitting param- 
eter for the simulation of the converse magnetoelectric 
effect in the following. 



VI. EXPERIMENT VERSUS SIMULATION 

By calculating the magnetization values M c and M a 
according to Eq. (2) the converse magnetoelectric effects 
in the Fe 3 4 (38nm)/BTO hybrid (cf. Fig. 3) can be 
simulated for different magnetic field strengths assum- 
ing that ferroelastic domain reorientation is the domi- 
nant contribution [cf. Eq. (1)]. The normalized magne- 
tization projection along the y direction as a function of 
the applied electric field obtained in this way is shown 
in Fig. 8, together with the experimental results taken 
from Fig. 3. An excellent agreement between experi- 
ment and simulation could be obtained by assuming a 
magnetoelastic coupling efficiency of 55% (x = 0.55). 
The fact that this value is well below 100% indicates 
that the magnetoelastic behavior of the Fes04 thin film 
is quite different from the bulk. Nevertheless, as obvi- 



ous from Fig. 8, the variation of the magnetization with 
the applied electric field can be well simulated for mag- 
netic fields jhqH > 50 mT using the two-region model de- 
scribed in Eq. (1). At magnetic field values close to the 
coercive field, converse piezoelectric effects and magnetic 
domain effects come into play, leading to large changes 
of the magnetic coercivity as a function of the electric 
field. 44 ' 46 These changes can not be explained within the 
simple two-region model alone. In this case, the two- 
region model must be further extended to take into ac- 
count piezoelectric effects as well as strain effects due to 
polarization rotations to simulate the electrical manipu- 
lation of the magnetic coercive fields. 



VII. CONCLUSION 

We have investigated strain-mediated converse mag- 
netoelectric effects in extrinsic multiferroic hybrid struc- 
tures consisting of epitaxially grown ferrimagnetic Fes04 
thin films on ferroelectric BTO substrates at room tem- 
perature. As evident from the time evolution of the in- 
tensity recorded by in-situ RHEED, the growth of Fe304 
thin films onto BTO substrates is highly disordered for 
the first few monolayers and then continuously evolves 
into a true layer- by-layer growth mode. SQUID magne- 
tometry measurements on these Fes04 thin films reveal 
a saturation magnetization close to the bulk value as well 
as a pronounced Verwey transition, which is a hallmark 
of the quality of Fea04. Furthermore, butterfly-shaped 



changes of the magnetization as a function of the applied 
electric field were observed at room temperature. The 
use of miscut BTO substrates, where only a single type 
of ferroelastic a domains exist, allows us to quantitatively 
simulate the observed converse magnetoelectric effect by 
using a simple two-region model. This model considers 
only effects caused by reorientations of ferroelastic c do- 
mains into ferroelastic a domains in the BTO substrate. 
In the framework of this model, the magnetization values 
of those parts of the Fes04 thin film elastically clamped 
onto ferroelastic c domains and a domains are considered 
as electric field independent quantities, which can be cal- 
culated using magnetoelastic theory. By comparing the 
results of the simulations to the experimental data, we 
demonstrate that the room-temperature converse mag- 
netoelectric effects of Fe304/BTO-based multiferroic hy- 
brid structures can be well understood on the basis of 
the two-region model for applied magnetic fields larger 
than the saturation field of Fes04. The detailed under- 
standing of the converse magnetoelectric effect at room 
temperature in BTO-based multiferroic hybrid structures 
opens the way to develop novel extrinsic multiferroic hy- 
brids exhibiting giant electric-field induced changes of the 
magnetization. 
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